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"o ■ Abstract 

The isobaric yield ratios of the fragments produced in the neutron-rich 48 Ca and 64 Ni projectile 



fragmentation are analyzed in the framework of a modified Fisher model. The correlations between 
the isobaric yield ratios (R) and the energy coefficients in the Weiszacker-Beth semiclassical mass 
formula (the symmetry-energy term a sym , the Coulomb-energy term a c , and the pairing-energy 
term a p ) and the difference between the chemical potential of neutron and proton (fi n — u p ) are 
investigated. Simple correlations between (u n — u p )/T, a c /T, a sym /T, and a p /T (where T is 
the temperature), and InR are obtained. It is suggested that (fi n — u p )/T, o c /T, a sym /T, and 



a p /T of neutron-rich nuclei can be extracted using isobaric yield ratios for heavy-ion collisions at 
intermediate energies. 
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I. INTRODUCTION 



The construction of a new generation of radioactive nuclear beam facilities has stimu- 
.ated much research into isospin physics Q. In heavy-ion reactions at intermediate energy, 
multifragmentation of the reaction system is generally observed in violent collisions and 
there is evidence that both subsaturated and supersaturated densities can be explored in 
such collisions |2j-|4|. Work in this area has concentrated on exploring the nuclear equation 
of state (EOS) and the liquid-gas phase transition in nuclear matter 5|-|8|. Isotopic yields 
in heavy-ion collisions provide a good probe for studying the nature of the disassembling 
nuclear systems. Many studies on fragment emission have attempted to use fragment yield 



distributions, either singly or by comparison to those of similar reactions, to explore the 



metry energy of the emitting source at different densities and temperatures |2|, |4j, [a, [7|, |9 



sym- 



20]. 



The nuclear symmetry energy of a finite nucleus is an important parameter in the EOS of 
an asymmetric nucleus and in various process in astrophysics and nuclear astrophysics. But 
the symmetry energy is difficult to measure experimentally and there are large differences 
in the theoretical results between different models and even within the same model with 
different parameters [l|. 

In a recent work analyzing isobaric yields [2l| , the ratio of the symmetry-energy coeffi- 
cient to temperature, a, s , ; m/T, as a function of fragment mass A was studied in a modified 
Fisher model (MFM) 22|, |23|. The Coulomb-energy coefficient to temperature (a c /T), and 
the pairing-energy to temperature (a p /T) were also extracted at the same time. For the 
symmetry-energy term, the extracted values from experiments are in good agreement with 
those calculated for the final fragments in the ground states. The pairing effect is clearly 
observed in experiments and strongly supports the hypothesis that the observed effect origi- 
nates at the end of the statistical cooling-down process of the excited fragments. A compar- 
ison between the Coulomb coefficients extracted experimentally and those calculated shows 
significant differences. 

In this article, on the basis of the theory of a modified Fisher model 22|, |23(, which 



was adopted in Ref . 21] , the correlation between the logarithm of the isobaric yield ratio 
lni? and a sym /T, a c /T, a p /T and (/i n — ft p )/T for fragments produced in 140 A MeV 48 Ca 
+ 9 Be and 64 Ni + 9 Be reactions (experimental data are taken from Refs. 24|, |25|) are 



analyzed. Coefficients of the volume energy, the surface energy, the Coulomb energy, the 



symmetry energy and the pairing energy in the Weiszacker-Beth semiclassical mass formula 
and (p n — p p )/T for fragments will be extracted using these correlations. 



II. ISOBARIC YIELD RATIOS IN THE MODIFIED FISHER MODEL 



23] , the yield of fragment with mass 



Following the modified Fisher model theory p, |22| . 
number A and I = N — Z, Y(A, I) is given by 

Y(A, I) = CA- T exp{[W(A, I) + p n N + p p Z]/T + Nln(N/A) + Zln(Z/A)}, (1) 

where C is a constant. The A~ T term originates from the entropy of the fragment, and 
the last two terms are from the entropy contributions for the mixing of two substances in 
the Fisher droplet model [26|. p n and p p are the neutron and proton chemical potentials, 
respectively, and W{A, I) is the free energy of the cluster at temperature T. W(A, I) is 
given by the generalized Weiszacker-Beth semiclassical mass formula 27|, |28| at a given 
temperature T and density p: 

W(A, I) = -E sym - a c (p, T)Z(Z - 1)/A^ 3 + a v (p, T)A - a s (p, T)A 2 ' Z - S(N - Z), (2) 

where the indices v, s, c, and sym represent volume, surface, Coulomb, and symmetry ener- 
gies, respectively. The symmetry energy E sym , can be divided into a volume-symmetry 

' p,T ) + S s (p,T), in which 



sym- 



s, 



term S v and a surface-symmetry term S s , i.e. E. 

S v (p,T) = a vsym (p,T)I 2 /A and S s (p,T) = a ssym (p,T)I 2 A 2 / 3 b 



pirical mass formulas, the pairing energy S(N, Z) is given by 31] 



301 ] . Following the semiem- 



S(N,Z) 



(3) 



' a p {p,T)/A 1/2 (odd-odd), 

(even-odd), 

— a p {p ) T)/A 1 ' 2 (even-even), 

The yield ratio for fragments, R(I + 2, /, A), between isobars differing by 2 units in I is 
defined as 

R(I + 2, 1, A) = Y{A, I + 2)/Y(A, I) 

= exp{[^(J + 2, A) - W(I, A) + (p n - p p )}/T + S mvx {I + 2, A) - S mix (I, A)}, (4) 

where S m i X (I,A) = Nln(N/A) + Zln(Z/A). To simplify the description, the density and 
temperature dependence of the coefficients in Eq. (J2J) is written as o, = aj(p, T) (where 



% = v , s, c, ssym, vsym and p represent the volume energy, the surface energy, the Coulomb 
energy, the surface-symmetry energy, the volume symmetry energy, and the paring energy, 



respectively). The temperature dependence of a^ at low T has been studied [29]. 
Inserting Eq. (J2j) into Eq. flU), one gets 

R(I + 2, I, A) =exp{[(^ n - n P ) - 4a sym (I + I) /A + 2a c (Z - I) /A 1 ' 3 

-S(N + 1,Z-1) + S(N, Z)]/T + A(J + 2, I, A)}, (5) 

where A(J + 2, 1, A) = S m i x (I + 2, A) — S m i X (I, A). Similarly, one can define the fragment 
yield ratio R(I + 4, /, A) between isobars differing by 4 units in / following Eq. (J3J) as 

R(I + 4, 1, A) = Y(A, I + 4)/Y{A, I) 

= exp{[W(I + 4, A) - W(I, A) + 2(fi n - n P )]/T + S rmx {I + 4, A) - S mix (I, A)}, (6) 

and inserting Eq. (J2J) into Eq. (jSJ), one gets 

R(I + 4, J, A) =exp{[2(/i n - // p ) - 8a s?/m (/ + 2)/A + 2a c (2Z - 3)/AV3 

- 6(N + 2, Z - 2) + <J(iV, Z)]/r + A(J + 4, /, A)}. (7) 

Equations ([6]) and ((Tj) assume that a s and a^ are the same for isobars, and omit the surface- 
symmetry-energy term as Ref. [2l|. In this case, a vsyrn is written as a sym according to Ref. 



2l|. 



For isobars with I = — 1 and 1 = 1 (which are mirror nuclei), A(l, —1,A) = 0, and the 
contributions from the symmetry term and the mixing entropy term in Eq. (jSJ) drop out 
and, the pairing term also cancels out because these isobars are even-odd nuclei. Taking the 
logarithm of the resultant equation, one obtains 

ln[i?(l, -1, A)] = [(/i n - fi p ) + 2a c (Z - l)/A l / 3 ]/T. (8) 

Following Eq. (j4]), an isobar with odd "I" is an odd-even nucleus and the pairing energy 
is zero, so one gets 

ln[R(I + 2, 1, A)] = [(/i„ - n P ) - 8a sym /A + 2a c {Z - lj/A^/T + A(I + 2, 1, A). (9) 

Considering ratios of isobars with "1-2", "I", and "1+2", and assuming isobars with "I" in 
R(I, I — 2, A) and isobars with 1+2 in R(I + 2, 1, A) are isotopes, one can reach 

(8a sym /A+2a c /A 1/3 )/T = \n[R(1, 1-2, A)] -ln[i?(/+2, 1, A)]-A(1, 1-2, A) + A(I+2, 1, A). 

(10) 



In Eq. ([IP]) , isobars with "1-2" are taken as the reference nuclei. 

Following Eq. ([6]) and taking isobars with "I" as the reference nuclei, one gets 

ln[i?(I +4,1, A)} = [^ n -fx p )-8(I + 2)a sym /A + 2a c (2Z-3)/A^ 3 ]/T + A(I + A,I,A). (11) 

Taking isobars with "1-2" as the reference nuclei, the difference between 21n[i?(/+2, J, A)] 
and hi[R(I + 4, /, A)} can be written as 

(8a sym /A+2a c /A 1/3 )/T = \n[R(I+2, 1, A)} -ln[R(I+A, /, A)} - A(J+2, /, A) + A(/+4, J, A). 

(12) 
It can be found that 

\n[R(I + 2, /, A)] - \n[R(I + 4, /, A)} + \n[R(I + 4,1 + 2, A)} 

-A(I + 4, / + 2, A) - A(I + 2, 1, A) + A(J + 4, 1, A) = 0. (13) 

Similarly, taking isobars with "1-2" as the reference nuclei, the difference between 
2\n[R(1, 1-2, A)] - \n[R(I + 2, 1, A)} reads 

l(fin-V P )+2a c Z/A 1/3 ]/T = 2\n[R(1, 1-2, A)]-\n[R(I+2, 1, A)]-2A(1, 1-2, A)+A(1, 1+2, A). 

(14) 



In Ref. [21], [(/x n — /z p )]/T for different reaction systems is expressed as [(/i n — fj, p )/T]i = 
[(^n — Hp)/T]o + Afi(Z/A)/T, in which (Z/A) = {Z v + Z t )/(A P + A t ), and p and t represent 
the projectile and the target nuclei. Taking one reaction system as the reference and fitting 
the different reaction systems, Afi(Z/A)/T for each reaction system can be fixed. Then 
o-sym/T and a c /T are extracted from the values of (fi n — [i p )/T for each reaction system. 
Here, a sym /T, a c /T and (/i n — ^ P )/T can be extracted using Eq. ffTQj) . ffT2|) . and ( JT4|) only 
for one projectile fragmentation reaction. 

Using Eq. ((H]), ( ITUl) . ( fl2|) . and f fT41) . we analyze the yield ratios of isobars produced 



in the 140A MeV 48 Ca + 9 Be and 64 Ni + 9 Be reactions |24j, [25|. In Fig. [T] and EJ the 
correlations between \n[R(1, 1 — 2, A)] — \n[R(I + 2, 1, A)} and the mass number (A) of 
fragments [Eq. ffTQj) ]. and the correlations between ln[R(I + 2, 1, A)] — ln[i?(J + 4, /, A)] and 
A of fragments [Eq. f)12p ] are plotted. These correlations are fitted using a function y = 
(8a sym /A + 2a c /A l ^)/T, in which A is the argument and a sym /T and a c /T are parameters. 
The values for ln[i?(l, — 1, A)] — ln[i?(3, 1, A)] shows a large difference to the values for 
\n[R(I,I - 2, A)] - \n[R(I + 2, 1, A)}. For isobars with / > 3, the values for \n[R(I,I - 
2, A)} - ln[R(I + 2, 1, A)} overlap. 
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FIG. 1. (Color online) Correlations between (8a sym /A + 2a c /A 1 ' 3 )/T for different isobars with / 
and A of fragments in the 140 A MeV 48 Ca + 9 Be reaction. The lines are the fitting results using 
Eq. (HDD. 
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FIG. 2. (Color online) Correlations between (8a sym /A + 2a c /A 1 ' 3 )/T for different isobars with / 
and A of fragments in the 140 A MeV 64 Ni + 9 Be reaction. The lines are the fitting results using 
Eq. (HO]). 



In Fig. El and @J the correlations between 2\n[R(I, 1-2, A)} - \n[R(I + 2, I, A)} and 
2Z/A 1 / 3 of isobars (using Eq. [T4l are depicted. The correlations are fitted using a function 
y = [([An — jjp) + 2a c Z/A 1 / 3 ]/T. The values for [(/i n — /x p ) + 2a c Z/A 1 / 3 ]/T for isobars with 
"I=-l" and "/ =1" shows large differences, but these decrease and there is an overlap as / 
increases. 



In Fig. [51 the correlations between ln[i?(l, — 1, A)] — ln[i?(J + 2, J, A)] for isobars with / 
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FIG. 3. (Color online) Correlations between 2ln[R(I, 1-2, A)] - ln[R(I + 2, /, A)] and 2Z/A 1 / 3 of 

fragments produced in the 140 A MeV 48 Ca + 9 Be reaction. The lines are the fitting results using 

Eq. flU). 
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FIG. 4. (Color online) Correlations between 21n[fl(J, 1-2, A)] - ln[R(I + 2, 1, A)} and 2Z/A 1 ? 3 of 
fragments produced in the 140 A MeV 64 Ni + 9 Be reaction. The lines are the fitting results using 
Eq. d). 



and 2Z/A 1 / 3 (using Eq. (TT5l) ) for fragments are displayed. The correlations are fitted using 
Eq. dHJ. The values for [(fi n - /x p ) + 2a c Z/A 1 / 3 ]/T for isobars with "I = -1" and "I = 1" 
show large differences, but these decreases and there is an overlap as I increases. 

Omitting the difference between (fj, n — /j, p )/T for nuclei with different /, and taking isobars 
with "I = -1" as the reference nuclei, the difference between ln[i?(l, — 1, A)]—ln[R(I+2, 1, A)] 



can be written as 

(8Ia asym /A + 2a c /A 1 / 3 )/T = hx[R(l, -1, A)} - ln[R(I + 2, /, A)] + A(J, / + 2, A). (15) 

In Fig.0 the correlations between ln[i?(l, — 1, A)] — ln[i?(J+2, J, A)] for isobars with different 
I and A for fragments are plotted. The correlations are fitted using Eq. (fl~5l) . The values for 
ln[i?(l, — 1, A)] — ln[R(I + 2, /, A)} for different isobars increase as the I of isobars increases. 
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FIG. 5. (Color online) Correlations between ln[i?(l, — 1,A)] — ln[R(I + 2, 1, A)] for isobars with 
different I and A of fragments in the 140 A MeV 48 Ca + 9 Be reactions. The lines are the fitting 
results using Eq. ([T5]) . 



Following the same methods in Ref. 2l|, the pairing term a p /T for isobars with 1=0 and 
1=2 are rewritten here: 

a p /T ~ (sgn)-A^ 2 {hi[R{2, 0, A)]--{]n[R(l, -1, A)]+ln[R{3, 1, A)]-A(3, 1, A)}-A(2, 0, A)}, 



(16) 



and for isobars with 1 = 2 and 1 = 4, 



a p /T ~ ( S ^)V/ 2 {ln[i?(4, 2, A)]-^{ln[i2(l, -1, A)]-31n[i*(3, 1, A)]+3A(3, 1, A)}-A(4, 2, A)}. 

(17) 
Similarly, one can have 

a p /T ~ ( S ^)V/ 2 {ln[ J R(2, 0, A)]-i{ln[J2(l, -1, A)]-31n[i2(3, 1, A)]+A(3, 1, A)}-A(2, 0, A)} 

(18) 



and for isobars with 1 = 2 and I = 4, 

a p /T ~ ( S ( ? n)iA 1 / 2 {ln[i?(4, 2, A)] - ^{ln[R(3, 1, A) - 31n[i?(5, 3, A)]] 



-A(3, 1, A) + 3A(5, 3, A)} - A(4, 2, A)}, (19) 



and 



a p /T ~ ( S ^)^ 1/2 {hi[i?(4, 2, A)] - ~{ln[i2(3, 1, A)] + ln[i2(5, 3, A)] 

+A(3, 1, A) + A(5, 3, A)} - A(4, 2, A)}. (20) 

For an (odd, odd) nucleus sgn = 1 and for an (even,even) nucleus sgn = — 1. The approxi- 
mations assumed in Eq. (fl6|) and (fT7|) are (/i n — /i p )/T, a sym /T and a c /T in ln[i?(3, 1, A)] and 
ln[i?(4, 2, A)] and are the same as those in ln[i?(l, — 1, A)]. Similar approximations are made 
in Eq. (ITSj) . (lT9~j) and (I20J1 . In Fig. El correlations between a p /T and A of fragments produced 
in the 140 A MeV 48 Ca + 9 Be reactions are plotted. a p /T of isobars with 1 = and 1 = 2 
are extracted using Eq. f fT6|) -f l20|) . respectively. In Eq.(ll6|). the chemical term, symmetry 
term and Coulomb term of R(2, 0, A) are assumed to be equal to those of R(l, —1, A). In 
Eq. ( IT8|) . (/i n — /j, p )/T, a sym /T and a c /T of R(2,0,A) are assumed to be equal to those of 
R(3,1,A). For (even,even) isobars with 1=0 and 1=2, the extracted a p /T using Eq. ( TT6l) 
is bigger than that using Eq. (IT5|) . while for (odd, odd) isobars the extracted a p /T using 
Eq. ffi~6]) are smaller than those using Eq. ( ITS]) . 

Due to lack of data for cross sections of mirror nuclei in the 64 Ni projectile fragmentation, 
figures like Figs. [5] and |6] are not plotted for 64 Ni. From Figs. [Uto|6j it can be seen that the 
correlations between isobaric yield ratios and (/i n — /j, p )/T,a syrn /T, a c /T, and a p /T can fit 
the measured data well. But for isobars with big J, there are not enough data to form chains 
and it is impossible to extract the values for a sym /T and a c /T for these very neutron-rich 
isobars. 

The method discussed above has the great advantage that the analysis can be performed 



in a single reaction and there is no need to calibrate Afi(Z/A) as in Ref. |2l|. The extracted 
values of (/x n — /x p )/T, a sym /T, a c /T, and a p /T are at a specific temperature associated with 
the incident energy. (/i n — /j, p )/T,a sym /T, a c /T, and a p /T are all temperature dependent. 
To study the dependence of (//„ — [i p )/T, a sym /T, a c /T, and a p /T on temperature, projectile 
fragmentation at different energies should be investigated. 




FIG. 6. (Color online) Correlations between a p /T and A for fragments produced in the 140 A MeV 
48 Ca + 9 Be reaction, (a), (b), (c), (d), and (e) are results obtained using Eqs. (fl~6|) . (fT7|) . ([18]) . 
flUD, and ([UD, respectively. 

III. SUMMARY 

In summary, the coefficients of the Coulomb energy a c /T, symmetry energy a sym /T, 
pairing energy a p /T, and (/i n — [i p )/T have been studied by analyzing the yield ratios (R) 
of isobars in projectile fragmentation in the framework of the modified Fisher model. Very 
simple correlations between (/i n — fi p )/T, a c /T, a sym /T, a p /T, and R are obtained. It is 
found that these correlations can fit the experimental results well and can be used to extract 
the symmetry energy of the neutron-rich nuclei. 
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